Mean excitation numbers due to anti-rotating term (MENDART) in cavity QED 

under Lindbladian dephasing 
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We study the photon generation from arbitrary initial state in cavity QED due to the combined 
action of the anti-rotating term present in the Rabi Hamiltonian and Lindblad-type dephasing. We 
obtain a simple set of differential equations describing this process and deduce useful formulae for 
the moments of the photon number operator, demonstrating analytically that the average photon 
number increases linearly with time in the asymptotic limit. 
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In the 2008-th paper by Werlang et al. fl'| a puzzling 
quantum effect was noticed from numerical simulations: 
when a two- level atom interacts with a single mode of the 
radiation field in a cavity by means of the Rabi Hamil- 
tonian, while subject to standard Markovian dephasing 
mechanism, the average intracavity photon number ex- 
hibits a linear growth with time. Such asymptotic photon 
generation due to decoherence occurs because for pure 
dephasing processes the environment may be viewed as 
a unmonitored detector making random nondemolition 
measurements of the number of quanta in the atom-field 
system [2, y|, whilst in [3-|y| it was shown that nonde- 
molition measurements can pump energy into the sys- 
tem via the destruction of quantum coherence provided 
the anti-rotating term is kept in the light-matter interac- 
tion Hamiltonian (i.e., without performing the Rotating 
Wave Approximation [3]). Besides, the pure dephasing 
reservoirs always possess a finite temperature (see, e.g. 
[2;, 3] for microscopic deduction) and hence store an infi- 
nite amount of energy, so the additional system energy is 
continuously supplied by the environment and the First 
Law of Thermodynamics is not violated (for the discus- 
sion concerning the Second Law of Thermodynamics in 
systems subject to frequent quantum measurements see 

Although the phenomenon of photon generation due to 
decoherence was explained qualitatively in [l|, |8| , no sat- 
isfactory analysis was carried out to derive analytically 
whether for the pure Markovian dephasing the average 
photon number de facto increases linearly with time and 
whether this growth saturates for large times. So the 
aim of this paper is to investigate analytically the behav- 
ior of Mean Excitation Numbers due to Anti-Rotating 
Term (MENDART), such as mean photon number and 
its variance or atomic excitation probability, and investi- 
gate their asymptotic characteristics in the simplest case 
of Markovian dephasing. We shall show that for any ini- 
tial state in the asymptotic limit the mean photon num- 
ber {n) indeed increases linearly with time, the average 
value of the photon number second moment (n^) grows 
quadratically with time, and the atomic excitation prob- 
ability Pe attains a constant value. So this paper pro- 
vides the missing mathematical explanation for the phe- 
nomenon of steady photon generation due to Lindblad- 



type decoherence in the presence of the anti-rotating 
term. 

Our starting point is the Markovian master equation 
for the density matrix p that takes into account both 
the atomic and cavity field phase-damping (dephasing) 

iiai 

P = -i[H, p] + -^[(JzP<7z - p) + lc {2npn - n^ p - pr?) , 

(1) 
where 7^ (7c) is the atomic (cavity) dephasing rate and 
H is the Rabi Hamiltonian ii3j[ll| (we set h=\) 

H = Lon+-az+9{a + a^){(y+ + (y-) (2) 

that includes the anti-rotating term {aa^ + a^ a^). Here 
a and a^ are the cavity annihilation and creation oper- 
ators, n = a^a is the photon number operator, and cj, 
ri and g are the cavity frequency, the atomic transition 
frequency and the atom-field coupling constant, respec- 
tively. The Pauli operators are defined as cr_ — \g){e\, 
c^+ = |e)(5l and a^ = |e)(e| - \g){g\, so that kets \g) 
and |e) can be interpreted as atomic ground and excited 
states, respectively. 

Expanding the density matrix in the Fock basis as 

00 
p = ^ (a„,„|g,n)(g,m| -t-5„,m|e,ri)(e,m| 

n,?n— 

+c„,m|5,n)(e,m| -hC„|e,n)(5,m|) , (3) 

where a„„i, 6„_m and Cn.m are time-dependent coeffi- 
cients, we obtain the exact set coupled differential equa- 
tions (the prime stands for the time derivative) 

ajj„ = i[uj{m-n)+i"fc{n~m) ]an^,n + ig{\/mcn,m-i 



- V^Cn^l + Vto + lc„,m+l - Vn+lcj;j_„+i)(4) 

- ^/nCn^l,rn + V^ + '^C*n+l,n - V" + lc„+l,ml5) 
^n.m — ^Jn,m^n,m ' ^dyw ^^^n^ra — 1 \y) 



where fn,m = w {m — n) + ft + {[ja + Jci'iT' — m)^]. In the 
strong dephasing Umit, (7a + 7c) > |g|, one expects on 
physical ground that due to the decoherence the terms 
Cn,m rapidly attain some constant values, so assuming 



that c' 



we get 



9 



■{Vn+lb, 



'n+l,r?i 



^^^n,m — 1 



In 

+^/r^fe„_l,„l - Vm + la„,,„+i) . (7) 

Now we substitute the expression for Cn,m back into 
equations (|ll)-(l5])[l3| and define new coefficients a„,m — 
e-*"*(™-")a„,™ and 6„,^ = e'^'^'^'^-^'fon,^ that are 
slowly varying functions of time. Assuming that \g\ <S 
Ci;,n (a condition that holds in cavity QED experiments 
unless the so-called 'ultra-strong coupling regime' |12l.[l3j 
is achieved) we neglect the rapidly oscillating terms and 
obtain the following effective differential equations for the 
diagonal probability amplitudes 

a'n = -[{vi+V2)n + V2]a„ + [vinbn-i+V2in+l)bn+i] (8) 



K = -[{vi+V2)n+vi]bn + [v2nan-i+vi{n+l)an+i], (9) 
where d„ = a„^rn b„ = b„^n and we defined coefficients 
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275^ 



275' 



(Lj-nY + -f^ 



{lu + ny + 72 



(10) 



with 7 = 7c -|- 7a standing for the total dephasing rate. 

One can easily verify that the normalization condition 
is maintained, J2^=o(^n + Ki) = 0^ so the equations (|8])- 
(|9]) are consistent and lead to the following coupled dif- 
ferential equations for the low-order MENDART 

(n(i))' = V2 + {vi - V2) [Pe{t) + {n<J,{t))] (11) 



Pe{t)' ^V2- (t^i + V2) [Pe{t) + {na,{t))] (12) 

{na^t))' = V2 - 2 (vi - V2) {n{t)) (13) 

- (vi + V2) [Peit) + {na.it)) + 2 {n^a,{t))] 



{nHt))' - 2(c.2 + r)2+^2) 



(14) 



X [-fg''{l + A{n{t)))-uoU{na,{t))'\ 



These equations cannot be solved analytically due to the 
coupling to the dynamical variable {Ti?az{t)j which obeys 
another differential equation. 

However, one can deduce the general formula for the 
average photon number {n(t)) by noticing the similarity 
in the last terms of equations pT|) and (|12l) . One gets 



{n{t))-{niO)) 



^2 + fi2 _|_ ^2 



{g^-ft - ujn[P,it) - P,{0)]} 
(15) 



where Peit) < 1 is still an unknown function of time|20|. 
Furthermore, in the asymptotic regime i — >■ 00 we expect 
from the equation (fT2|) that Pe attains a constant value. 
Imposing limi_j.oo Peit)' = we obtain from P^ -([T ^ 



lini [Peit] 



(ncT.m 



ujfl 



uj^ + Q'^ + 72 



lim (na^it))' =0 



(16) 



(17) 



lim (n^a.it)) = -- 



2ujn 



Cj2 + ri2 ^_ ^2 1^^ 
and from equation (I14[) we get 



lim {n{t)) (18) 



lim 



(n^t))' 



273' 



uj^ + n^ + 72 



1 



4 lim 



it)) ■ (19) 



Hence, in the asymptotic regime i — > 00 we have the 
following rules for the Asymptotic MENDART (AMEN- 
DART, as coined in Q) for any initial state: a) (n(i)) 
and — (n'^az{t)\ increase linearly with time; b) On?{t)\ 
increases quadratically with time; c) Pe{t) and {ncTz{t)) 
attain constant values. We solved numerically the effec- 
tive differential equations ([S])-® and verified that the 
formula (jlSp is correct for all times, thereby accounting 
for the linear growth of {n{t)) noticed in [1] from numer- 
ical data, while the equations (fT6|) - (fT9|) agree with the 
numerical results in the asymptotic regime. 

In the figures [T]-[3] we compare the exact dynamics re- 
sultant from the original differential equations (HI) -([6]) to 
the effective dynamics governed by the simplified equa- 
tions ([S])-® for the parameters cj = 1, g = 4 x 10~2^ 
7a = 2g and 7c = 0. In the figure [T] ([2]) we consider the 
resonant regime SI = a; (dispersive regime Vi = lo — 2Qg) 
for the initial zero-excitation state \g, 0). We plot the dy- 
namical behavior of observables (n), (n.^), (n?(Jz), (ncTz) 
and Pe calculated from the original differential equations 
(I!])-®. Within the thickness of the lines these curves are 
indistinguishable from the graphs resulting out of the ef- 
fective equations ([S])-®. To exemplify the difference be- 
tween the original and effective differential equations, we 
show the zoom for the behavior of (n^) at initial times: 
the solid line depicts the exact dynamics, and the dashed 
line - the effective one. The observed discrepancies are 
quite small and appear because c^ ^ does not becomes 
zero instantly as was assumed in our analysis; neverthe- 
less, these minor differences are irrelevant regarding the 
asymptotic behavior. In the figures we also show the pho- 
ton number distributions calculated numerically at the 
time gt — 300 according to the original differential equa- 
tions (bars) and the effective ones (dots). Once again, 
the agreement is excellent. 

In the figure [3] we consider the initial state p(0) — 

Ptherm ® |e)(e|, whcrc ptherm IS the thermal state of the 

' Electromagnetic field whose photon number distribution 

is pn = n^ / {n -\- 1)""*"^, where n is the average photon 



number. We set Q — uj, n — 0.3 and show the asymptotic 
behavior of (n), (n^) and {r?u^ obtained from the orig- 
inal differential equations (figure |3K) and the zoom of (n) 
and Pe for initial times (figure [5}d, where the dashed lines 
represent the solutions of the effective differential equa- 
tions). We see that asymptotically the behavior agrees 
with equations P^ ^ dT^l) . although the transient dynam- 
ics cannot be reliably described by the equations dS])-®. 




FIG. 1: Exact and effective dynamical behavior of principal 
observables and the photon statistics for the time gt — 300 in 
the resonant regime, il = lj. 




FIG. 2: Same as figure[T]in the dispersive regime, Q — lj — 20g. 

Regarding the practical observation of the asymptotic 
linear photon growth inside the cavity due to decoher- 
ence, it seems quite unlikely in current cavity or circuit 



(a) 



(b) 




FIG. 3: Behavior of (n). 



^>, (^ 



and Pe in (a) the 



asymptotic regime and (b) during the transient regime for 
small times. The initial state is ptherm ® ie)(e| with average 
photon number {n(0)) — 0.3. 



QED implementations because the photons (and atomic 
excitations) would be lost due to radiative and nonra- 
diative relaxation processes. As example, let us consider 
the current state of the art circuit QED implementations. 
The typical parameters are [13] : ui ^ fl ^ 8 GHz and 
g ^ 0.3 GHz, while the dephasing rate is of the order 
of 7a ^ 1 MHz, although it can be made large at will 
(usually one desires to decrease 7 and not to increase 
it). Considering a high value for the total dephasing rate 
7 ^ IGHz the resulting asymptotic photon growth rate 
due to dephasing would be ~ 1 MHz. This value is of the 
same order of magnitude as the cavity relaxation rate for 
a rather high cavity quality factor Q ^ 10*, so (n(i)) 
would saturate at some (small) value instead of showing 
an asymptotic growth, as calculated explicitly in 8, 9, 16[ 
for standard quantum optical master equation. Some 
photons escape to the outside world via radiative dissipa- 
tion channel so they could be ultimately detected outside 
the cavity, but in this case different models predict dif- 
ferent photon emission rates depending on assumptions 
made about the reservoir [13, 13 (in particular whether 
it is Markovian or not). 

Recently a more sophisticated microscopic model was 
developed for deducing the master equation in the pres- 
ence of the anti-rotating term, valid in a specific regime 
of parameters [iJl- According to that model, the phe- 
nomenon of dephasing-induced generation of photons is 
greatly exaggerated by the Lindblad-type master equa- 
tion ([T}, and instead of the linear asymptotic growth 
the average photon number saturates at some value that 
strongly depends on the reservoir spectral density \13^. 
Nevertheless, the very phenomenon of photon generation 
due to decoherence persists and our formulae provide the 
upper bound for the photon generation rate. From the 



qualitative viewpoint, in realistic (lossy) cavity QED ar- 
chitectures this phenomenon would lead to a parameter- 
dependent heating of the system slightly above the ther- 
mal equilibrium values [9|, |l6| , depending on the atom- 
field detuning, coupling strength and the dephasing rate. 
In summary, we obtained simplified differential equa- 
tions describing the process of photon generation (from 
vacuum or any other state) due to the combined ac- 
tion of the anti-rotating term and the standard Lindbla- 
dian dephasing in Markovian cavity QED, whose validity 
was confirmed by extensive numerical simulations. From 
these equations we deduced analytical formulae describ- 
ing the overall behavior of MEND ART for arbitrary ini- 



tial state, demonstrating that asymptotically the mean 
photon number (n) increases linearly with time at the 
rate 2jg^{uj'^ + Q'^ + -f^)~^, (n^) grows quadratically with 
time and the atomic excitation probability attains a con- 
stant value. 
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Notice that the obtained asymptotic photon generation 
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although in that paper the mathematical approach was 
oversimplified. 



